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Gn . Abstract 
On 

^ \ A Langevin equation with multiplicative noise is an equation schematically of 

Q ' the form dq/dt = — F(q) +e(q) £, where e(q) £ is Gaussian white noise whose 

amplitude e(q) depends on q itself. I show how to convert such equations 
into path integrals. The definition of the path integral depends crucially on 
the convention used for discretizing time, and I specifically derive the correct 
[ path integral when the convention used is the natural, time-symmetric one 

that time derivatives are (qt — qt_At)/At and coordinates are (qj + q 4 _At)/2. 
^vq [This is the convention that permits standard manipulations of calculus on the 

action, like naive integration by parts.] It has sometimes been assumed in the 
literature that a Stratanovich Langevin equation can be quickly converted to 
a path integral by treating time as continuous but using the rule 9(t=0) = ^. 
O I show that this prescription fails when the amplitude e(q) is q-dependent. 

Oh! 
CD . 

43 ■ I. INTRODUCTION 



On 



Let £ be Gaussian white noise, which I'll normalize as 

(&(*)£;(*')> = ^(t-f')- 

It's long been known that a Langevin equation of the form 

^fc = -*i(q)+& (1.2) 
can be alternatively described in terms of a path integral of the formj] 

rq(t")=q" 



P(q",q',t"-t')= f q [dq(t)\ exp - / L(q,q 

Jq(i')=q' Jtf 



;i.s) 



1 For a review of background material in notation close to that I use here, see, for example, chapter 4 of 
ref. 0| . The most substantial difference in notation is that my F is that reference's -|f. 



Here, V(q", q', t) is the probability density that the system will end up at q" at time t if it 
started at q' at time zero. However, the exact form of L depends on the convention used in 
discretizing time when defining the path integral. With a symmetric discretization,^ 



P(q",qV)=jimiv/ a 
with Lagrangian 



q(t)=q" 
q(0)=q' 



L t 



exp 



qt - qt-At q< + qt-At 
At ' 2 



J (1-4) 



^(q,q) = ^|q + F| 2 -iF M . (1.5) 



Here and throughout, I adopt the notation that indices after a comma represent derivatives: 
Fij = dFi/dqj and Fijk = d 2 Fi/dqjdqk- N is the usual overall normalization of the path 
integral, which I will not bother being explicit about. 

What hasn't been properly discussed, to my knowledge, is how to correctly form such a 
symmetrically-discretized path integral for the case of Langevin equations with multiplicative 
noise (meaning noise whose amplitude e(q) depends on q). Schematically, 

d 

— q { = --Fi(q) + e ia (q) £ a , (1.6) 

with £ as before (|1 . 1| ) . I will assume that the matrix Ci a is invertible. There are a wide variety 
of applications of such equations, but I'll just mention one particular example of interest 
to me, which motivated this work and for which a path integral formulation is particularly 
convenient: the calculation of the rate of electroweak baryon number violation in the early 
universe 0. 

By itself, the continuum equation ( |1.6| ) suffers a well-known ambiguity. To define the 
problem more clearly, we must discretize time and take At — > at the end of the day. 
Specifically, I will interpret ( |1.6|) in Stratanovich convention, writing 



qt - fh-M = -At F(q) + e(q) (1.7) 
(Za£w)=&tn5a b 6 te , (1.8) 



where 

_ _ q< + q t _At 



;i.9) 



[The fact that I've labelled the noise £ t instead of £ t _At in (|1-7|) is just an inessential choice 
of convention.] The Stratanovich equation ( p..7|) is equivalent to the Ito equation 

q< -qt-At = -AtFi(q t _ A t) + e(q t -At) € t , (1.10) 



For a discussion of what changes if other discretizations are used in this case, see refs. 
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with 



- 

~ 77 e ia,j£ja- 



1.11) 



I will give two different methods for deriving the corresponding path integral. The result 



is 



„ , /-q(i)=q" 
' Af^O 7q(0)=q' 



x exp 



n deti 



qt + qt-At 



1 -1 



qt - gt-At qt + qt-At 
At ' 2 



;i.i2) 



£(q, q) 



2Q 



(q + F) t g l3 (q + F), 



1 1 -i • fi 



where 



fl'ij — ( e )ia( e )ja- 



;i.i3) 



;i.i4) 



This differs from a result previously given by Zinn- Justing [|TJ by the inclusion of the last 
term in L. Zinn- Justin's derivation was done in continuous time, resolving ambiguities using 
the prescription 9(t=0) = |, which is known to work in the case where e(q) is constant. 

Because of the confusion surrounding these issues, I will show how to derive the result in 
two different ways. First, I will follow the standard procedure for directly turning Langevin 
equations into path integrals, but I will be careful to keep time discrete throughout the 
derivation. The second method will be to start from the Fokker-Planck equation equivalent 
to the Langevin equations ( |1.7| ) and ( |1.10[ ) and to then turn that Fokker-Planck equation 
into a path integral, again using standard methods. 



II. DIRECT DERIVATION FROM THE LANGEVIN EQUATION 



Rewrite the discretized Langevin equation (|1.7| ) as 

E r - C = 0, 

where r is a discrete time index and 

E r = e-^q) [q T - q T _i + AtF(q)], 



(2.1) 



(2.2) 



with q ee (q T + q T _!)/2. The corresponding path integral is obtained by implementing these 
equations, for each value of r, as 5-functions, with appropriate Jacobian, integrated over the 
Gaussian noise distribution: 



ri(t)=i" 
P(q",q',t) = lim N \ 

At^O 7q(0)=q' 



Y[d£ T exp | 



2VtAt 



dq T S(E T - £ T ) 



det 



dE T ' a 



r'a;T"i \ dq T ' 



7 (2.3) 



3 Specifically, eq. (4.79) or ref. @. See also ref. §) and section 17.8 or ref. for a continuum time discussion 
of formulating the path integralior this problem using ghosts. 
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The noise integral then gives 



pww,*) - J [n^] -p(-^xx) r det, • 



In our case (|2.2|) , the determinant takes the form 





oqi 

oqi dq2 

dq 2 dq3 



dq3 dq4 



det 



det 



V 



J|det 



4 
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(2.4) 



(2.5) 



The registration of the diagonals is determined by the nature of the initial boundary condi- 
tion, which is that go is fixed. From fl2.2|) , we then have 



det ^ 

T'a;r"i \ C/g r //.j 



n dct 



% + 5 (e-^afc.i (?r ~ ffr-l)* + At |( e ~ lF )<V. 



J] det [e 1 (q r )j \ < ]Jdet + \ e ja (e (g T - g T _i) fe + At ±e ia (e 1 F) a J I , (2.6) 



where all e's and F's should now be understood as evaluated at q r . Now rewrite the 
determinants in the last factor of ( |2.6| ) as exponentials in the usual way, using 



det(l + A) = e tvln( - 1+A) = exp tr [l + A - \A 2 + 



(2.7) 



To construct a path integral, we need to keep track of the terms in each time step up to 
and including O(At), but we can ignore corrections that are higher-order in At. For this 
purpose, the size of q r — q r -i should be treated as 0(\/At), which is the size for which the 
g 2 term in the action [the exponent in (|1.12 )1 becomes 0(1) per degree of freedom. So, using 
the expansion (|2.7|), we get 



det 



Sij + \ e ja (e l ) aKi Aq k + At ~e ja {e 1 F) CI;i 



expj (constant) + \ e ia (e l ) ak>i Aq k + At |e io (e 1 F) Qii 

-| e ia (e- 1 ) afeim e m6 (e- 1 ) W)i Ag A A qi + 0[(At) 3 / 2 ] 



(2.8) 



It is the AgAg term in this equation, which came from the second-order term in the ex- 
pansion (|2.7|) , that will generate the difference with the result quoted in ref. HJ. Putting 
everything together, we get the path integral (|1.12| ) with Lagrangian 



Lm- q) = (q + F) igi j{q + F), - ^ F M + ^ e^e ka , k {q + F)< 



1 

2 



At 



+ — e ia (e 1 



ak,m&mb\e )bl,iQk<ll 



(2.9) 
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We can simplify this by realizing that the %qi in the last term can be replaced by its leading- 
order behavior in At. Specifically, recall that q t — qt_At is order y/At. So one can go for a 
large number of discrete time steps 1 <C iV <C 1/ At without any net change in q at leading 
order in At. Moreover, the force F doesn't have any net effect, at leading order in At, over 
that number of steps. The upshot is that %qi can be replaced at leading order in At by its 
average over a large number of steps, ignoring F and treating the background value of e(q) 
as constant. The Gaussian integral for q in ( |1.12| ) and ( |2.9| ) then gives the replacement rule 



Ml ^^(9'% (2.10) 
at leading order in At. This substitution turns ( |2.9| ) into the result (|1.13| ) presented earlier. 



III. DERIVATION FROM FOKKER-PLANCK EQUATION 



The Stratanovich Langevin equation ( |1.7| ) is well-known to be equivalent to the Fokker- 
Planck equation 



d_ 

dqi 



Q d 



d qj 



(3.1) 



where P = P(q, t) is the probability distribution of the system as a function of time. This 
is just a Euclidean Schrodinger equation, and one can transform Schrodinger equations into 
path integrals by standard methods. Specifically, rewrite the equation as 



with the Hamiltonian 



P = -HP, 



H = ^ pi e ia (q) pj e ja (q) -zp-F(q). 



(3.2) 



(3-3) 



To obtain a path integral with symmetric time discretization, it is be convenient to 
rewrite H in terms of Weyl-ordered operators. The Weyl order corresponding to a classical 
expression (9(p,q, t) is defined as the operator O w with 



(q\6 w \q'}= f e^'-^0(p,i(q + q'),t 



(3.4) 



For the sake of completeness, I'll briefly review how to obtain Weyl ordering of operators in 
simple cases by considering the application of the operators to an arbitrary function ip(q). 
For example, 
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Doing the p integrals with our Hamiltonian (|3.11|) then reproduces the path integral ( |1.4| 
with Lagrangian 



£(q, q) = ^ (<li + + l^e ia e k a,k) g i:j (qj + Fj + fOe^e^) + u(q). 



(3.15) 



Now note that 



(9 



(3.16) 



and so 



„ \&ia&ka,k)9ij\&jb&lb,l) ~\~ U &ia,i&ja,j T 



4 



1 

— 2 *>' "g" e * a >j e j < v' 



(3.17) 



Combining ( |3.15|) and ( |3.17|) reproduces the Lagrangian ( |1.13|) asserted in the introduction. 
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